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Abstract 



We show that the boundary between quark-dominated and hadron-dominated re- 
gions of nuclear structure may be blurred by multi-nucleon quark clusters arising from 
color percolation. Recent experiments supporting partial percolation in cold nuclei 
and full percolation in hot/dense nuclear matter include: deep inelastic lepton- nucleus 
scattering, relativistic heavy-ion collisions and the binding energy in 5 He\. 



1 Introduction 



Heavy ion experiments and lattice-gauge simulations of Quantum Chromodynamics (QCD) 
intensively investigate the quark-gluon phase transition relevant to the early universe. An- 
other challenge is to understand, under varying external conditions, the roles of fundamental 
quarks and gluons in hadronic and nuclear matter. Where are they essential to replace 
low-energy descriptions of nuclei based on hadrons (neutrons, protons and mesons)? From 
lattice QCD simulations, we understand quantitatively the masses of the hadrons. Similarly, 
we know that quark and gluon color charges make the strong coupling constant grow at dis- 
tances larger than the hadron size and generate their own confinement within color-neutral 
hadrons. However, a theoretical understanding of longer range QCD dynamics, beyond 
the size of a hadron, is still a challenge for lattice simulations. We present evidence for 
an intermediate-range phenomena, the formation of multi-quark clusters in cold and hot 
equilibrated baryonic matter that blurs the boundary between the hadron and quark-gluon 
phases. 

We adopt classical percolation theory (TJ |2j [3j 0] to bridge the gap between quark and 
hadron physics. In our hybrid approach we take nuclear matter to be composed of color 
neutral clusters, which consist at zero temperature of 3, 6, 9, 12, ... quarks and at finite 
temperature we include quark- ant iquark (qq) states for the mesons. These larger units have 
small probabilities in cold nuclei, but are conducting color, like subdomains of ionized plasma 
in a electrically neutral gas. We assume there is a common size parameter for the lowest 
mass nucleon and pion which determines the overlap of hadronic states at finite density and 
consequently the probabilities to form these multi-hadron clusters. We expect the presently 
neglected gluon content of these systems to become important at very short distances. 

In fig.l we summarize the percolation approach to nuclear physics where, between the 
well-established hierarchy of quark/gluon physics and nucleon/meson physics, we insert a 
new layer consisting of color neutral multi-quark clusters. Percolation theory allows us to 
connect these regions and appears sufficient to unify the microscopic and macroscopic physics 
with an accuracy of ~ 10%. Solutions for cluster probabilities in nuclei have previously been 
reported [HE]. Current experiments probe nuclei with high resolution and/or high excitation 
energies where this quark-clustering layer becomes manifest as we now discuss. Future 
electron-nucleus (e-A) experiments at JLAB with 12 GeV could enhance the resolution by 
an order of magnitude and reveal precise cluster details. 

1 Deep Inelastic Scattering on Nuclei with x > 1 

For deep inelastic lepton-nucleus scattering, in Born approximation, the Bjorken- variable 
x = 2 ® nU (Q is the magnitude of the 4-momentum transfer to the target, v is the laboratory 
energy loss of the lepton, and is the nucleon mass) measures the fraction of the struck 
quark's momentum relative to the nucleon momentum. Therefore, a constituent of an iso- 
lated nucleon carries x < 1. Quark clustering in nuclei leads to a simple prediction for the 
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Figure 1: (Color online) The separation of quark physics, the physics of color neutral multi- 
quark clusters and and nucleon/meson physics. 

region x > 1 where cooperative effects are required to provide scattering support. 

Following our previous work [3], we assume that the nucleus consists of 3, 6, 9, 12, 
... quark clusters with definite internal structure (i.e. defined structure functions) and 
cluster probabilities are governed by intermediate range correlations evaluated with realistic 
correlated wavefunctions for A < 4 nuclei, then scaled according to average density for 
heavier nuclei [I]. We omit effects of additional strong density fluctuations due to, for 
example, alpha-cluster formation in A > 4 nuclei. Two nucleons form a 6-quark cluster 
when they are separated by less than a critical distance d c = 2R C = lfm, where R c is a 
nucleon's critical "color percolation" radius. A third nucleon with separation less than d c 
from the previous two then joins to make a 9-quark cluster, etc. The probability Pi(A) that 
a quark selected at random in nucleus A originates from an i-quark cluster (for i — 6,9, 12...) 
decreases rapidly with increasing quark number. It follows that the e-A structure functions 
are determined by the 6-quark cluster for 1 < x < 2 and by the 9-quark cluster for 2 < x < 3. 
Independent of the cluster structure functions, per nucleon ratios of deep inelastic cross 
sections in two nuclei A and B will be constant in the respective x- intervals and reflect only 
the ratio of cluster probabilities of these nuclei, i.e. for 1 < x < 2 

(l/A)da/dxdQ 2 (eA) _ p 6 (A) 

{l/B)da/dxdQ 2 {eB) ~ p 6 {B) ' ^ 

and for 2 < x < 3 

(l/A)da/dxdQ 2 {eA) _ p 9 (A) 
(l/B)da/dxdQ 2 (eB) ~ p 9 (B) ' 

and so forth for larger clusters. We have shown how this model successfully describes the 
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inclusive inelastic cross sections on deuterium in the region 1 < x < 2 [6] when the 6-quark 
cluster component is included. 

The intense JLAB electron beam has allowed these suppressed regions of phase space to 
be measured pj|. In Fig. [2] we show the measured ratio (discrete points) [7] of per nucleon 
responses of 56 Fe to 12 C together with our predictions (horizontal bars) [4] in the 6-quark 
cluster region 1 < x < 2 and the 9-quark cluster region 2 < x < 3. We estimate the 
theoretical ratios have 10% uncertainty. In the former kinematic region the experimental 
ratio for 56 Fe to 12 C is 1.17 ±0.04 ±0.11 [7] in good agreement with our theoretical prediction 
of 1.17 ± 0.12. The experimental ratio above x = 2 is quotedjT] as 1.44 with about 33% 
uncertainty which agrees with our theoretical ratio of 1.39 ± 0.14. 
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Figure 2: (Color online) JLAB data for the ratio of per nucleon responses of 56 Fe to 12 C 
over a range of Bjorken xb = x that exceeds unity [7J, the limit for data on an isolated 
nucleon. We note there is a data point, r = 2.2 at x — 2.8, which is far off the vertical 
scale. Multi-quark clusters support a nuclear response above x — 1. The horizontal lines 
show the theoretical predictions (1.17, 1.39) for regions dominated by (6,9)-quark clusters 
respectively [1]. 

Although we derive encouragement from the agreement (within large uncertainties) for 
the ratio in the 9-quark cluster region we believe this agreement may be accidental for 
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two reasons. First, the large x values in these experiments come with rather small energy 
transfer v in the rest system of the nucleus. If one translates these v values into a time 
resolution then one finds that this "snapshot" has a rather long exposure time, namely 
about At = l/0.3GeV = 0.7 fm/c. Therefore, the short-lived fluctuation of a 9-quark cluster 
may be represented only partially in the 2 < x < 3 cross section. Future e-A experiments 
with new or upgraded facilities can elucidate the dependence on exposure times. They 
would also allow higher momentum transfer Q, which suppresses the coherent effects of 
quasi-elastic quark-cluster knockout [3] at x = 1,2, etc., which may be visible in Fig. [2j 
Second, the calculations [I] were performed without three-nucleon (NNN) interactions since 
they were unavailable at that time. NNN correlations should reflect the repulsive NNN force 
and suppress 9-quark probabilities overall. Due to their intermediate range, we expect the 
NNN force will also lead to non-trivial A-dependent effects. These complications motivate 
significant theoretical effort to modernize the quark cluster probability calculations. 

An alternative interpretation [S] of these deep inelastic scattering data employs short- 
range two-nucleon (NN) correlations to generate very high lGeV/c) nucleon relative 
momenta so that the struck quark may obtain a large momentum and exceed the kinematical 
boundary x = 1. The high momenta are obtained by strongly repulsive short range NN- 
potentials acting at separation distances of r < 0.7 — 0.8/m. In so far as this nucleon-based 
approach generates similar kinematics, it may appear to represent a non-relativistic picture 
dual to our quark-cluster approach. However, color-conducting multi-quark clusters offer 
both a relativistic approach and additional experimental consequences. For example, virtual 
photons at high resolution can peer into the quark cluster zone of color conductivity and 
measure details of the quark cluster structure functions. We may discover, for example, that 
percolating quarks carry less of the total momentum sum rule than isolated nucleons even 
though x > 1 regions are populated. This seems plausible on the basis of elementary physics 
since the color confining region is enlarged in the 6-quark cluster [9]. 

A recent experiment probing the x < 1 region suggests that, while 3 He, 4 He and 12 C 
responses follow a density scaling trend, the 9 Be response does notpjl]. In addition, a detailed 
study of hadron-nucleus data at x > 1 shows significant deviations from simple density 
scaling for 6 Li and 7 Li nucleipT]. We understand these Li and Be anomalies as motivating 
the inclusion of alpha-clustering effects which are expected to be especially significant for 
these nuclei. For this reason and our neglect of NNN potentials, we prefer to focus on larger 
systems that we believe are less sensitive to these corrections. 

2 Relativistic Heavy Ion Collisions 

We now show that the quark cluster description is appropriate at high excitation energies 
produced in relativistic heavy ion collisions. Cluster probabilities are mainly a function 
of the local hadron density and, to proceed to finite temperature, we introduce a pion 
component. We assume that the critical separation distance for the pions to transition 
into larger quark clusters is the same as for nucleon percolation, d c = lfm, and that d c is 
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temperature independent. Therefore we adjust the baryon chemical potential so that the 
sum of the baryon and pion densities is the same as the density of baryons at T = to 
maintain approximately constant color percolation probability. Then 

Phadrons Pbaryons Ppions 

(p„T). (3) 

We set the pion chemical potential fi n = at T = to fix the reference baryon chemical 
potential /i& at normal nuclear matter density 0.17 baryons//m 3 (Fermi momentum kf = 
1.36 fm~ l ). The nucleon dispersion relation is influenced by the nuclear mean field potential 
for which we take [12] 

e(k) = a/ (k 2 + (m N + V s (k)f) + V {k) (4) 
where the scalar and vector self-energies are momentum dependent 

V s (k) = -0.35GeV(l -0.05(k/k f ) 2 ) (5) 
V (k) = +0.29GeV(l-0m(k/k f ) 2 ) (6) 

The freeze-out data [13] extracted from heavy ion collisions at high temperatures are 
shown in Fig. 3 as a function of T and pb along with the theoretical (dotted) curve which 
corresponds to a mostly hadronic system of pions and nucleons. On this curve, each quark 
is percolating 22% of the time in i > 3-quark clusters. We also plot in Fig. 3 a curve with a 
higher percolation rate which we call "percolating nuclear matter" and find pb where more 
than 50% of all nucleons have transitioned to 6, 9, 12, ... -quark clusters at T = 0. This T = 
state requires pbaryon = 0.46/m -3 i.e. 2.7 times nuclear density and therefore has a larger 
chemical potential than nuclear matter. For this case, we also adopted p n = 120MeV, 
used to fit the transverse momentum dependence of the pion spectra from a kinetically 
equilibrated pion distribution [13]. Taking into account also the A- resonance at 1.236 GeV 
shifts the higher temperature percolation curve (solid line) towards the curve corresponding 
to nuclear matter density (dotted line). 

The percolation model is a geometrical model and is not sensitive to the decay properties 
of resonances, but more to the number of hadronic states in a certain volume. The classical 
percolation transition is a cross over transition in finite systems and the transition itself is 
not well defined. Thus quark percolation blurs the boundary between the hadron phase and 
the quark/gluon phase. 



3 Extended Quark Orbitals in Hypernuclei 

Another avenue to explore quark cluster structures is to replace one of the nucleons with 
a strange spectator, a A— hyperon, producing a hypernucleus. The successful explanation 
of the hypernucleus spin-orbit potential based on quark dynamics [15J motivates a closer 
look at possible quark cluster effects. In particular, we focus on a longstanding puzzle [16] 
of the sudden decrease of the binding energy in 5 HeA compared with the binding energy 
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Figure 3: (Color online) The boundary between hadron- and quark/gluon-matter from our 
quark percolation model for 50% percolation rate (upper curves) and 22% percolation rate 
(lower curve) in a temperature T versus baryon chemical potential phase diagram. The 
parameters of the three theory curves are specified in the legend. The arrow indicates \ib for 
nuclear matter at T = 0. The freeze-out data are compiled from different experiments [13]. 

of a baryon-based calculation with two-body forces. The 4 He nucleus is doubly magic in 
the language of the shell model with protons and neutrons but the full occupation of the 
0s -state by protons and neutrons does not preclude the addition of a A— hyperon to the 0s 
-state. With two-body forces alone one then obtains overbinding of the A=5 hypernucleus by 
1 — 3MeV. One way out is through effective many-baryon potentials. Indeed, calculations 
with purely baryonic degrees of freedom and a three-body A-NN force resulting from an 
explicit E- component in the wave function [T7] seem to show the repulsive character for the 
total three-body contributions. However, if the three quark substructure of the A matters, 
i.e. up(u)+down(d)+strange(s), then for the 12-quark cluster configuration of 4 He, the A's 
up and down quarks cannot occupy the Os-shell which can accommodate only 12 u- and d- 
quarks according to the Pauli principle. The observation of such extended quark orbitals 
would be the best proof for color percolation in nuclei as G. Baym has already indicated in 
1979 pp. 
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Let us make a simple estimate of this effect. The mean spacing in the constituent quark 
shell model would be larger by a factor three based on the ratio of nucleon to constituent 
quark masses 
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Including a size factor for the smaller 12 quark cluster compared to the nuclear size could 
easily boost this mean spacing by a factor of two. 
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Figure 4: (Color online) Sketch of the two possible views of the Lambda attached to the 
the lowest configuration of two neutrons(n) and two protons(p). On the left, the Lambda 
and the nucleons are structureless baryons with spin projections indicated. On the right, 
their quark substructure is taken into account with spin and color attributes sketched. In 
the latter case, the 12-quark cluster configuration that appears with probability pu = 0.007 
in the quark cluster model, forces the up(u) and down(d) quarks of the Lambda to be in 
the next higher, Pauli-allowed, orbit 0p q . The Lambda's strange(s) quark may occupy the 
lowest 0s q level. 

Therefore a 12-quark cluster with two extra u- and d- quarks would need an extra 
excitation energy of approximately 150-300 MeV. Our quark cluster model provides jl] 
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Pu^He) = 0.007, a small probability for such an extended 12-quark cluster in 4 He (Fig. 4). 
Multiplying this probability with the Pauli-required u- and d- quark excitations of 150-300 
MeV would give a decrease in 5 He^ binding of 1 — 2MeV, in rough agreement with experi- 
ment. This one case is not sufficient to establish the existence of extended quark shell model 
states. One should look at heavier hypernuclei where the additional A has to split up in a low 
lying s-quark and an excited ud-pair. Future precision experiments could explore the binding 
energy systematices of heavier closed shell systems such as: 16 0a, 17 0a and 16 N^, 17 N\. 

Additional experimental tests of quark clustering can be performed such as measuring 
high mass di-lepton pair production from nuclear targets [T8] (Drell-Yan process). 

If one asks whether our understanding of QCD is complete, the question may be recast 
to whether additional theoretical or experimental work can improve the situation decisively. 
Compressed baryonic matter [19] is an area where our theoretical and experimental under- 
standing is not as well developed as baryon-free matter. For nuclei a theory which contains 
nucleons and multi-quark cluster states is an efficient way to describe the transition from 
purely nucleonic matter to quark matter at high baryon density. A hybrid model with quark 
clusters would define a major advance if one would be able to rewrite the smooth infrared 
limit of lattice QCD in a field theoretic continuum picture. Such a picture should also con- 
tinue the phase diagram into the totally percolated phase where perhaps the formation of 
a qq -dimer condensate [20] and BCS-pairing [21 J of quarks can be found in analogy to the 
physics of cold atoms. Quark percolation and these other features of the QCD phase diagram 
would also impact astrophysics. 
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